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Abstract. We present the results of a search for physics beyond the Standard Model (BSM) using data of 
1.1 fir 1 integrated luminosity collected by the CMS experiment at the LHC. Fully hadronic final states were 
selected based on the "stransverse" mass variable M T 2 and interpreted in various models of supersymmetry 
(SUSY). Two complementary analyses were performed targeting different areas of the SUSY phase space. All 
backgrounds were estimated using both simulation and data-driven methods. As no excess of events over the 
expected background was observed exclusion limits were derived. 



1 Introduction 

We describe a search [ 1 1 for physics beyond the Standard 
Model in pp collisions collected by the Compact Muon 
Solenoid (CMS) detector [2] at the Large Hadron Collider 
(LHC) at a centre-of-mass energy of 7 TeV. The results are 
based on a data sample of 1.1 ftr 1 of integrated luminosity 
collected in 2011. We use the "stransverse mass" variable 
Mj2 [3| to select new physics candidates out of fully ha- 
dronic events. We divide our search into two channels: one 
targets high squark and gluino masses with a high Mji cut, 
the other heavy squarks and light gluinos with a medium 
M T 2 cut but including a b-tag and high jet multiplicities. 
In the following we describe the properties of Mji in sec. 
[2] our analyis strategy and event selection in sec. [3] and 
the background estimation in sec. [4] In sec. [5] we state the 
results of our search, and draw a conclusion in sec. [6] 



2 The search variable M T2 

The variable Mji was introduced to measure the mass of 
primary pair-produced particles where both particles decay 
into detected and undetected particles (e.g. the lightest su- 
persymmetric particle (LSP)). It is a generalization of the 
transverse mass mj in case of two identical decay chains 
each containing unobserved particles. The variable Mji is 
defined as 

Mj2(yny) — min [max ( m.9}, m? )1 , (1) 

where mj is the transverse mass of the visible system and 
the corresponding LSP x of the decaying sparticle: 

mf = ^( m ^(0)2 + m 2 + 2 (Ef (i) Ef° - P ; ,s(,) • p* (0 ). 

In this analysis the stransverse mass Mji is not used for 
mass measurements but rather as a discovery variable (4). 

In order to associate all visible decay products to the 
decay chains of the two sparticles we cluster the jets of an 
event into two "pseudojets" using a hemisphere algorithm 
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described in [5|, Sect. 13.4. As seeds (inital axes) the di- 
rection of the two (massless) jets are chosen which have 
the largest invariant mass. We then associate a jet k to the 
pseudojet i rather j if the Lund distance is minimal: 

Ei Ej 

(Ei — PiCosOii) ~ < (E : — p icosOik) ~. 

v F, < k \ E . + Ek) 2 y J 'J Jk '(Ej + E k ) 2 

2.1 Advantages of M T2 

In order to gain a better understanding of the behaviour of 
Mj2 we take the case where we set all masses to zero and 
assume no initial-state radiation (ISR) or upstream trans- 
verse momentum 1 . In this simple case M T o becomes 

(M T2 f=2pf m pf {2 \\ + cos<Pn), (2) 

where p"j il) is the transverse momentum of pseudojet i, 
and 0i2 the angle between the two pseudojets in the trans- 
verse plane. It can be observed that for symmetric events 
(p™ = pj'^) with large acoplanarity Mj2 behaves like 
the missing transverse momentum (MET). Thus SUSY with 
expected large MET will accumulate in the high M T 2 re- 
gion. However, back-to-back systems or balanced events 
will populate the region with small Mji- Thus Mji is ro- 
bust against QCD jet mismeasurements: Mismeasurements 
along one of the pseudojets results in Mti ~ GeV, while 
for asymmetric mismeasurements still Mri < MET. 



3 Analysis Strategy and Event Selection 

In this analysis we have established two search channels in 
order to be sensitive to different regions in the SUSY phase 
space. One approach, the High Mji analysis, targets events 
resulting from heavy sparticle production which is charac- 
terized by large MET and Mji- The second approach, the 
Low Mji analysis, is designed to be sensitive to the region 
where squarks are heavy and gluinos relatively light. Here 

1 Upstream transverse momentum is the transverse momentum 
which is not clustered to the pseudojets (e.g. jets outside of ac- 
ceptance). 
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gluino-gluino production is dominant, the gluinos giving 
rise to three-body decays with small MET. Also, as stops 
and sbottoms are expected to be relatively light, these events 
can be enriched with b-quarks. Thus the two strategies re- 
quire two different sets of selection cuts stated in table [T] 

Table 1. Event selection cuts which are specific to their strategies. 



High M T2 


Low M T2 


at least 3 jets 2 

H T > 600 GeV 3 
M T2 > 400 GeV 


at least 4 jets 
at least 1 b-tag 
H T > 650 GeV 
M T2 > 150 GeV 



Besides this channel specific selection we also require: 

- Lepton (<?, fi) veto to reduce W+Jets and it background. 

- min A0(MET, any jet) > 0.3 to reduce further QCD. 

- |MHT - MET | < 70 GeV to minimize the influence of 
unclustered energy (e.g. ISR) to the M T2 shape 4 . 

- MET tail cleaning cuts (e.g. noise filters) to filter out 
events with unphysical MET. 

For the selection of data we require the data to pass Mj 
trigger paths. 

4 Background Estimation Strategy 

For each type of background data-driven estimation meth- 
ods have been designed: QCD is estimated from the bulk 
of the Mt2 distribution as described in sec. |4.1| In order 
to reduce the effect of signal contamination and statistical 
fluctuations the electroweak and top background is esti- 
mated from an adjacent control region in Mj 2 . The pre- 
diction is taken from data in the control region scaled by 
Monte-Carlo (MC) ratio of the event yield in the signal 
region over the yield in the control region. Similarly the 
uncertainties are scaled by a MC ratio. The control region 
for High M T2 is defined as 200 GeV < M T2 < 400 GeV, 
for Low Mj-2 it is defined as 100 GeV < M T2 < 150 GeV. 

4.1 QCD background estimation 

The QCD estimation method is based on the two variables 
Mj2 and A0 min = minA0(MET, any jet). These two vari- 
ables are strongly correlated, but a factorization method 
can be applied if the functional form for the ratio r(M T2 ) - 
N(A(p min > 03)/N(A<f> min < 0.2) is known. From simula- 
tion studies it has been found that for Mj 2 > 50 GeV the 
ratio falls exponantially: 

N(Ad>- > 0.3) 
N(A<f> min < 0.2) 

This behaviour was confirmed in data. The estimate has 
been performed by a fit from data in the QCD dominated 
region of 50 GeV < Mj 2 < 80 GeV to extract the parame- 
ters a and b. In order to get also parameter c from data its 
value was fixed to the value of the ratio at M T2 = 200 GeV 
where the ratio still falls exponentially. 

2 The jet selection requires p T > 20 GeV, \rj\ < 2.4 

3 H T is the scalar sum of all jet-p r . 

4 MHT is the negative vectorial sum of all jet-p^. 



4.2 Z — > w background estimation 

Z — > vv is an irreducible background because the produced 
neutrinos leave the detector unmeasured and thus generate 
real MET. The number of Z — > vv+Jets events passing the 
event selection can be estimated from W — > yt/v+Jets via 

N Zvv (est) = W(+iv) R zw (4) 

^acc^recoj iso 

where W(fiv) is the number of W — > fiv events passing the 
event selection with additionally requiring one muon, Rzw 
is the ratio of Z — > vv events to W — > fiv events, e acc is the 
acceptance, and e reco /iso is the combined reconstruction and 
isolation efficiency. In order to reduce the it background 
in the W — > fiv selection a b-tag veto has been applied, 
while the residual it background has been estimated from 
the ^-tagged region. Furthermore e reco /i SO has been calcu- 
lated from Tag & Probe studies on Z — > U events while the 
acceptance and the ratio Rzw are taken from MC. 

4.3 W and Top background estimation 

The background due to W and Top has two sources: Ei- 
ther a lepton (e, fi) from a W has been unobserved due 
to acceptance cuts, or has been "lost" due to failing ei- 
ther identification, isolation, or reconstruction criteria. The 
other source are W decays into neutrinos and taus which 
decay hadronically. 

The number of events with a "lost" lepton has been 
estimated from the number of events with one lepton found 
in data. This number is then corrected for the probability 
to loose a lepton via the formula 

N pass veto = ^eco _ rfftlZ** , ( 5 ) 

where N r e e ™ is the number of events containing a lepton, 

N^n is the expected background from processes other than 
W or Top, and s e ^ is the probability for a W — > lv(l = e,fi, 
or t — > e,fi) passing all selection and reconstruction cuts. 

The number of events with hadronic tau decays are 
taken from simulation and validated by data: The W —> Iv 
(I = e, fi, t) kinematics in simulation has been validated in 
data with one muon. Furthermore, it has been shown that 
tau decays are well modelled in the simulation [6| justify- 
ing the use of simulation. 



5 Results 

The M T2 distributions for the High M T2 analysis and Low 
Mj 2 analysis are shown in Figs. [T] and [2] table [2] summa- 
rizes the results of the two analysis strategies. As no excess 
over background has been found limits have been set. 

5.1 Exclusion Limits 

First, model independent limits oncrx BR within our ac- 
ceptance has been derived by computing a 95% upper limit 
on the number of events using a CL S formulation |7 1. These 
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High M T , Analysis CMS Preliminary.*^ = 7 TeV, L = 1.1 fb' 




Fig. 1. M T2 distribution for the High M T2 analysis. The MC back- 
ground is normalized to 1.1 fb -1 . A possible SUSY signal (LM6) 
is overlayed. Data are shown as dots on top of the background. 
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Fig. 3. Combined exclusion limit in the mSUGRA/CMSSM 
(mo, m 1/2 ) plane with tan/? = 10. 



Low M .. Analysis CMS Preliminary.Vs = 7 TeV, L = 1.1 fb' 1 
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Fig. 2. M T 2 distribution for the Low M T2 selection. The MC back- 
ground is normalized to 1. 1 ftT 1 . A possible SUSY signal (LM9) 
is overlayed. Data are shown as dots on top of the background. 
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Fig. 4. Model pp — » gg, g — > bbx° for the Low M T2 selection as a 
function of the mass parameters m LS p and m s : Signal efficiencies 
(left), 95% CL upper limit on cross section of the model (right). 



6 Conclusion 



Table 2. Expected yield in the signal region from Standard Model 
(SM) background in simulation and from data-driven background 
predictions, as well as data yield for both analysis strategies. 





SM-MC 


Background prediction 


Data 


High M T2 


10.6 


12.6 ± 1.3 (stat) ±3.5 (syst) 


12 


Low M T2 


14.3 


10.6 ± 1.9 (stat) ±4.8 (syst) 


19 



Table 3. Observed and expected limits on cr x BR within the ac- 
ceptance ot the two analysis strategies. 





<tx BR (pb) 




observed limit expected limit 


High M T2 


0.010 0.011 


Low M T2 


0.020 0.014 



limits are shown in table [3] Exclusion limits at 95% C.L. 
have been determined in the mSUGRA/CMSSM (m , mi/ 2 ) 
plane. The results are shown in fig.[3]for An = 0, p > and 
tan/3 = 10 combining the High and Low MT2 selections 
by taking the best expected limit in each point. Besides 
the exclusion in the mSUGRA/CMSSM plane the results 
are interpreted in a so-called Simplified Models topology. 
This is a simple signal model with exactly one decay mode 
which is only constrained by the kinematics and the masses 
of the participating particles. The Low Mji analysis is in- 
terpreted in a model where gluinos are pair produced and 
each gluino decays into two fo-quarks and a neutralino, the 
LSR The limits on the model cross sections and the signal 
efficiencies are shown in fig.|4] 



We conducted a search for supersymmetry in hadronic fi- 
nal states using the Mji variable calculated from massless 
pseudojets. A data set containing 1.1 ftT 1 of integrated lu- 
minosity in yfs — 7 TeV pp collisions recorded by the 
CMS detector during the 201 1 LHC run was analyzed. Two 
complementary analyses were performed to probe a larger 
SUSY phase space. In both analyses the tail of the Mji is 
sensitive to a possible SUSY signal. As no evidence for a 
signal was found, we set upper limits on the cross section 
times branching ratio within our acceptance. Exclusions 
limits were established in the mSUGRA/CMSSM parame- 
ter space, as well as in a Simplified Model topology. 
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